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Abstract 

It is known that with restrictions on the type of the constitutive equations, 
Maxwell's equations in non-uniform media can sometimes be reduced to two 2nd or- 
der differential equations for 2 scalar quantities only. These results have previously 
been obtained in two quite different ways. Either by a "scalarization of the sources", 
where the relevant scalar quantities are essentially vector potential components and 
where the derivation was limited to isotropic media; or alternatively by using the 
"scalar Hertz potentials", and this method has been applied to more general media. 
In this paper it is shown that both methods are equivalent for gyrotropic media. We 
show that the scalarization can be obtained by a combination of transformations 
between electric and magnetic sources and gauge transformations. It is shown that 
the method based on the vector potential, which previously used a non-traditional 
definition of the vector potentials, can also be obtained using the traditional defini- 
tion provided a proper gauge condition is applied and this method is then extended 
from isotropic to gyrotropic media. It is shown that the 2 basic scalar Hertz po- 
tentials occurring in the second method are invariant under the source scalarization 
transformations of the first method and therefore are the natural potentials for ob- 
taining scalarization. Finally it is shown that both methods are also equivalent with 
a much older third method based on Hertz vectors. 

1 Introduction 

Usually problems in electromagnetics are reduced to solving a 2nd order vectorial equation 
for either the electric field or for the magnetic field [HE]- However for some problems, the 
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radiation of a dipole in a stratified medium being a good example, it is still useful to use a 
representation for the fields in terms of auxiliary functions (vector and scalar potentials or 
Hertz vector potentials) and then solving equations for these auxiliary functions instead 
of for the fields directly. It is well-known that for simple uniform electromagnetic media, 
Maxwell's equations can in this way be reduced to 2 scalar uncoupled 2nd order differential 
equations, corresponding with the TE/TM modes. For uniform media this has been 
generalized to more complex (decomposable) media [3]. Also for non-uniform, usually 
stratified, media such a scalar decomposition has been obtained. 

One of the earliest systematic treatments of this "reduction problem" was given by Nisbet 
[H E]. First a general formulation in terms of 2 Hertz vector potentials was given valid 
for non-uniform anisotropic media and then it was noticed that considerable freedom ex- 
ists for choosing these potentials since they can be subjected to a gauge transformation. 
Using this freedom the Hertz vectors could be reduced to single component vectors and 
conditions were derived under which the resulting differential equations for these compo- 
nents were of 2nd order, at least for isotropic media. The same idea was later extended 
to anisotropic media [SI [7]. The reduction to two scalar potentials was also extended to 
gyrotropic media [HI [9] and to even more complicated media [lO]. Only rather recently the 
case of a uniaxial medium was given explicitly [IT]. Whereas in these earlier publications 
[HI [7] a general coordinate system was considered, the later extensions to more complex 
media usually considered a cartesian coordinate system only and the non-uniformity was 
limited to a stratification along e.g. the symmetry axis of the uniaxial medium [TT]. Al- 
though the 2 scalar potentials are referred to as "scalar Hertz potentials" the link with the 
Hertz vector potentials is not obvious anymore and for the stratified uniaxial medium the 
(initially 4) scalar functions are instead defined by applying a Helmholtz-decomposition 
to the electric and magnetic field components perpendicular to the symmetry axis [TT] . 

Subsequently Weiglhofer and Georgieva [12], [13] arrived at the same scalar equations, at 
least for isotropic media, following a completely different method, which rests mainly on 
the so-called scalarization of sources. It was shown that arbitrary current density distri- 
butions can be replaced by equivalent distributions but oriented along a fixed direction. 
With a proper (unconventional) choice of the vector potentials the latter could then also 
be scalarized and the resulting equations are exactly the same as those found using the 
scalar Hertz potentials. This remarkable correspondence was noticed but no explanation 
was given [T3] . 

The main purpose of this paper is to shed some light on this finding, which can not be a 
coincidence. At the same time we will extend the source scalarization method explained 
in [111 [13] to more general (gyrotropic) anisotropic media. We believe that "source scalar- 
ization" can be understood best as an application of the well-known equivalence between 
electric and magnetic charges and currents. This is first presented in § [2 In the following 
sections the scalarization problem is solved using different potentials. As in most of the 
referenced papers we consider only cartesian coordinates. As medium a stratified initially 
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uniaxial medium is considered where the symmetry axis is perpendicular to the layers 
everywhere. The latter condition is necessary to avoid mixing between the longitudinal 
and transverse field components by applying the constitutive equations. In §[6] the theory 
is extended to a stratified gyrotropic medium. We will use c as a unit vector along the 
symmetry axis and c as the corresponding coordinate whereas transversal vector compo- 
nents will be labeled by ±, in particular the transversal nabla operator will be written as 



2 Source transformations 

We use Maxwell's equations in the standard form including electric and magnetic charge 
and current densities, where the latter are labeled by a superscript star 

— dB — * , , 

VxB = -^-J (1) 

— dD - 
VxH = — + J (2) 

V-D = p (3) 
V-B = p* (4) 

The possibly position dependent constitutive properties of the medium are given hy D = 
e ■ E and B = p ■ H. We will write the electric charge and current densities in general as 

J = ^ + xm (5) 
ot 

p = -V-p (6) 

where p and m are either given polarization and magnetization densities or must be 
considered as stream potentials for given p and J [5]. In either case p and m can be 
subjected to a gauge transformation which leaves p and J invariant [5] 

p' = p + V xG (7) 

m = m — — + Vg (8) 

where G, g are arbitrary functions. Since our aim is to transform transversal sources into 
longitudinal ones, 2 possibilities arise. With V ±g = —Tn±_ a transversal magnetization 
is turned into a longitudinal one gc. And with a longitudinal G = GcC a transversal 
polarization pj_ = — V^Gc x c is turned into a longitudinal magnetization — ^^c. 
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It is well-known that the polarization p and the magnetization m can equally well be repre- 
sented by magnetic charge and current densities, which is most easily seen by rearranging 
the terms in the Maxwell equations as follows 

Vx(E + e-^-p) = - — (5-/i-m)--^r — + V x e"^ ■ p (9) 

ot dt 

— d — 

Vx{H-m) = —{D + p) (10) 

V-(D+p) = (11) 
V-(B-fi-m) = -V-(/i-m) (12) 

The equivalent magnetic sources are thus given by 

p* = -V ■ (/X ■ m) (14) 

Contrary to the gauge transformation ((7]) (E]) , in this case the source transformation (from 
electrical charges to magnetic charges) is accompanied by the following field transforma- 
tions 

e' = E + e-^-p (15) 
b' = B-fi-m (16) 

In what follows we will label a polarization/magnetization density which is represented 
by magnetic charges by a superscript star {p*,Tn*). We can thus freely exchange p (or 
m) for p* (or m*) and vice versa as long as p + p*(or m + m*) remains invariant and as 
long as we take the field transformations f |T5ll(fT6l l into account. These "magnetic" stream 
potentials can also be subjected to a gauge transformation [5] 

e-'-p'* = + (17) 

p.fn'* = fi-m* -V xL (18) 

for arbitrary L, I. With I we can turn a transversal polarization into a longitudinal one and 
with L = LcC we can turn a transversal magnetization into a longitudinal polarization. 
For further reference the different representations are tabulated in [H 

Using the electric/magnetic charge transformations and gauge transformations if needed 
we can now scalarize an arbitrary current densit}ll] along a fixed direction defined by the 



^We will consider an electric current density, but the same method can be applied to a magnetic 
current density. 
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Table 1: The equivalent source contributions due to (external) polarization and magneti- 
zation. The starred quantities allow to make a distinction between sources modeled with 
electric charges/currents and those modeled with magnetic ones. 





polarization 


magnetization 


J 


tfp 
dt 


V X m 


p 


-V -p 




r 


-V X e"^ ■ p* 


dt 






—V ■ /i ■ m* 



unit vector c. The goal of the scalarization process is to replace the current density by 
equivalent electric and magnetic current densities parallel to c. For a stratified medium 
c is perpendicular to the layers and as mentioned in the "Introduction" this is also the 
direction of the symmetry axis of the uniaxial medium. The current density can always 
be written as 

J = j^c + V±v X c + V±u (19) 

In the context of scalar Hertz potentials the functions u, v are known as auxiliary functions 
[TI] and they can be found by solving the 2-dimensional potential problems 

V^M = V^-J± (20) 
Vlv = -c ■ (Vi X 7^) (21) 

Scalarization of the 2nd term in ( fT9l) is straighforward, since according to [H it can be 
attributed to a magnetization fn = vc which can also be represented by a magnetic 



current density along c 



— * da ■ m* dv , , 



where /i// is the permeability along c. This transformation is accompanied by a field 
transformation according to (fTBD 

B' = B-fi//vc (23) 



Scalarization of the last contribution in (1191 ) cannot be obtained simply by transforming 
electric into magnetic sources. However this transformation can always be combined with 
a gauge transformation. The last contribution in (fT9l) can then be scalarized by attributing 
the current density V_lm to a transversal polarization density p = V_l / udt which, as we 
have seen, can be turned into a longitudinal one in the "magnetic" domain by a proper 
choice of / in (fTTi ) (with L = 0) namely 

l = - J ^dt (24) 
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J'c = -e//- ( - 1 (25) 



where V_Le_L = has been used. We then end up with a scalarized current density 

d f u 
dc \e± 

Due to the 2 electric/magnetic transformations preceding and following the gauge trans- 
formation we must take into account a transformation of the electric field according to 

m 

e' -E = Vl = V (^j ^d?j (26) 

These results as well as similar ones for the magnetic current density are tabulated in [2] 
and [HI To conclude this section we make 2 remarks 

1. It can be proved that the scalarizations summarized in [2] and [3] are unique; 

2. Under these source transformations the total longitudinal current densities ( Jc + 
and J* + are invariant, since 

Jc + ^ = (VxF).c (27) 
J: + ^ = -{VxE).-c (28) 

and h' — H and E — E are either parallel to c or equal to a gradient. 



3 Vector potentials 

We introduce conventional vector and scalar potentials A, and also comparable poten- 
tials A , (f)* for handling the magnetic sources 

— — OA* 

B = V X A-fi- - /i ■ V0* (29) 

E = -^-V^-e-' -V X A (30) 
ot 

Substitution into Maxwell's curl-equations gives initially 

£(e,/x)A + e-V^ = 7 (31) 
C{fi,e)A* = T (32) 
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Table 2: Using proper transformations between electric and magnetic charge represen- 
tations from [T] and possibly gauge transformations of the stream potentials, arbitrary 
transverse electric current densities can be replaced by equivalent electric or magnetic 
current densities along the symmetry direction c. Each transformation is also accompa- 
nied by a transformation of the fields, shown in the last 2 rows. 



initial 
current 
density 


J = V±v X c 


J = V±u 


initial 
stream 
potenial 


fn = vc 


p = V± J udt 


final 
stream 
potential 


fn* = vc 




final 
current 
density 






e' -E 






b'-b 


-fl//VC 
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Table 3: Using proper transformations between electric and magnetic charge represen- 
tations from [T] and possibly gauge transformations of the stream potentials, arbitrary 
transverse magnetic current densities can be replaced by equivalent electric or magnetic 
current densities along the symmetry direction c. Each transformation is also accompa- 
nied by a transformation of the fields, shown in the last 2 rows. 



initial 
current 
density 


J* = V_Lf * X c 


f = V^M* 


initial 
stream 
potenial 


— ^ ^ — 

p = — e//f c 




final 
stream 
potential 


p = —e//v*c 


m* = -^( f ^dt 

ac \J /ix 


> 


final 
current 
density 


~f _ dv* -p. 


J* = -^^//^c 




E' -E 


V*C 




b'-b 
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where the operator £(e, /i) is defined by 

£(e,/i) = e- — + V x/i-1- Vx (33) 

The scalar potentials are eliminated using gauge conditions. If these gauge conditions 
do not mix-up the electric and magnetic quantities then also the resulting equations will 
remain uncoupled. We first mention the gauge conditions used by Nisbet [5] 

V-e-A + a^ = (34) 
ot 

V-fi-A* + a*^ = (35) 

where a, a* are scalars which can still be chosen. As will become clear further on these 
gauge conditions allow scalarization for a homogeneous medium only. This restriction is 
eliminated in [T2l [T3] but only for an isotropic medium by using a different definition for 
the potentials and using different gauge conditions. At first this leads to equations for 
A and A which are coupled but under the restrictions V^e = and V_L/i = these 
equations become uncoupled and scalarizable. Although this method can be extended 
to uniaxial media we prefer to stick with the conventional decompositions (l29l)(l30l) and 
we will now derive gauge conditions which allow to scalarize (1311) (l32l ) for a non-uniform 
uniaxial medium. 

We assume that the sources have already been scalarized so that in fl3TI)(!32ll only longitudi- 
nal current densities Jc, J* occur. Splitting these equations into longitudinal and transver- 
sal components the latter equations will only allow the null solution for the transversal 
components of the vector potentials if their longitudinal components do not occur in 
these transversal equations. These conditions are easily found by assuming A = Aji and 
equating the transversal component of the LHS of (ISTl) to zero 







_dc 





(36) 



and a similar "magnetic" equation, which is obtained by replacing unstarred quantities by 
starred ones and by switching the roles of e and /i. Obviously using the gauge condition 
(l3il) . this condition can only be met if e^, fi± do not depend on the longitudinal coordinate 
c and then only by choosing a = e//e±fi±. For a more general result we must instead 
choose the following gauge conditions 

V-fil^A + e±^ = (37) 

V-el^r+zix^ = (38) 
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From the longitudinal components of f l3Tl)(l32ll and these gauge conditions we then find 
the final scalarized equations 



>Cs(e,/i) 



■ Al* 



-J' 



-J' 



(39) 
(40) 



where scalarized quantities are now explicitly marked by an accent and with the (scalar) 
wave operator defined by 



d _ d d"^ 
£,(e, /x) = Vi + e//—el'— - e//fi^ — 



(41) 



Taking into account ([29l) (l30l) and the field transformations due to the prior source scalar- 
ization, the fields are given by 



H 
E 



dA'l 



dt 



c + /ij X c - V (fil^ J 



u 



c - eT V I A'! X c - V e 



u 



d_/AZ 
dc \ e_L 

dc V/i± 



dt 
dt 



(42) 
(43) 



4 Hertz vectors 



The main idea behind the use of Hertz vectors is to introduce an additional differentia- 
tion in such a way that the new potentials (= Hertz vectors) are governed by equations 
with the polarization/magnetization (or the stream potentials) as sources instead of the 
current densities. Since for such Hertz vectors the starred/unstarred stream potentials 
are equivalent we can eliminate e.g. the "magnetic charge" sources from the start and we 
can use instead of fl29l) (l30l) the simpler decompositions 



B + jj, ■ m 
E 



-1 — * 
e ■ p 



V X A 



(44) 
(45) 



still complying with ([T]) and (l4|) . The Hertz vector equations are obtained most easily in 
the temporal gauge [14j (0 = 0) and we then define the Hertz vectors Hg, H^ following [5] 



A 



' + ■ V X n. 



(46) 
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After substituting these equations in ([2]) ([3]) and following a standard procedure [E] we 
find the equations 

C{e,fi)ne = p + r + e-V^e (47) 

C{fi,e)Ilm = fi- (m+m*) + fl-Vtprn (48) 

where the operator >C(e, /i) has been defined in (l33l) . The functions 4'e,i'm can be chosen 
arbitrarely and are in fact redundant due to the gauge transformations (TTl) (l8l) (fT71) (fTSll . 
Substituting (Il6ll in (liHl and (Il5l) and replacing the 2nd order time derivative using (IT7I) 
and fl33l) we find the following symmetric expressions for the fields 

B + fi ■ rn* 

e-E + p 

Using the gauge transformations (TTl) (l8l) (fT71) (fTSll one can try to simplify the equations 
(ITfj ) (llSll ■ One possibility is to eliminate the magnetizations (fn + m* 0) so that also 
Hm 0. This is the essence of the 1-Hertz- vector method followed by Sein [H], in 
this case limited to a uniform isotropic medium and (therefore) using the conventional 
Lorentz gauge. Another possibility is to scalarize the stream functions so that p + p* —>■ 
p'ji and fn + fn* —>■ m'lc. For these scalarized sources it is now possible to choose the 
(gauge) functions ipe,ipm. in (IT7ll(ll8l) in such a way that ne,!!^ are also scalarized and 
thus have only components along c. To that end we assume Yie/m = ^{e/m)cC and collect 
the transversal components of e.g. ( l47l l 



QH _ 

V X — ^ + V X e"^ ■ V X n„ 

ot 



V X ■ V X He - V X 



dt 



(49) 
(50) 



^ [/il^ (V^nj] = exVxV'e (51) 
The transversal components of the vector potentials will then vanish if we choose 

V-e = — 7^ 52 

e_L oc fi± 

The longitudinal part of the same equation is given by 
and using (l52l l we obtain 

=-p'c (54) 
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where the scalar operator Cs has been defined in (|4T]) . A similar equation can be found 
starting from ( l48l l 

£,(/i,e) j j = -/i//m'* (55) 



Unlike the scalarized current densities, which are unique, the scalarized stream potentials 
and the corresponding Hertz vector components are not unique, since they can always be 
subjected to a gauge transformation. However we can in particular choose the scalarized 
stream potentials as follows 

Jc = ^ ^c = /^//^ (56) 

where J'c and J'* are the (unique) scalarized current densities. Comparing (I54l)f l55il with 
(I39l)fl40l) we conclude that in that case 

This correspondence can also be checked by comparing the field expressions (Il9l) (l50ll with 
(l42l ) f l43ll where in (l49l) (l50l l one should also take into account the field transformations due 
to the scalarization of the current densities (see [2] and [3]) . 



5 Scalar Hertz potentials 

The "scalar Hertz potential" formulation introduced by Weiglhofer [TT] starts by decom- 
posing the fields and the equations into transversal and longitudinal components and 
parts. From these equations the longitudinal components Ec, He can be eliminated leav- 
ing 4 equations for the 4 unknown transversal components. Up to this point the method 
is identical to the 4x4 matrix method used for solving Maxwell's equations in strati- 
fied media [17]. However for dealing with the source terms and unlike the 4x4 matrix 
method, the transversal field components are then expressed using scalar (Hertz) potential 
functions 

Ei_ = V±$ + V± X Gc (58) 

H± = v±n + v±xm (59) 

In what follows we will give a compact derivation of the scalar Hertz potential equations, 
following [11]. Splitting Maxwell's curl-equations parallel and perpendicular to c we obtain 

c-{V^xH^) = Jc+^ (60) 
c-(VxxE^) = -{J: + ^) (61) 
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and 



V.ff. - ^ = ^ X (7x + ^) (62) 
Taking the cross-product with c, the latter 2 equations become 

The 4x4 matrix method is based on (l63l) and (l64ll . without the source terms, and where 
(l60l) and (IBTl) are used for eliminating E^ and i^c, after inserting the constitutive equations 

D, = e//E, D^ = ei_E^ (66) 
B, = fi//H, B^ = fiJI^ (67) 

Weiglhofer [H] deals with the source terms by operating with V_l- on the 4 equations 
(l62l ) - ([65l) and by introducing the auxiliary functions already defined in (|20|) and (|2T1) . 
Using the constitutive equations (166]) ([671) and the decompositions ([58l)(l59l) all terms then 
contain the laplacian V^^ which can be dropped yielding 

On the other hand the longitudinal equations ( i60ll (IBTl) become 

-Viv^ = Jc + e//^ (72) 

Vie = + (73) 



Finally dEc/dt and dHc/dt can be calculated from (l68l) - ([7Tl) as a function of 6 only 
and when substituted in (r72l)(l73l) one obtains two uncoupled 2nd order equations in 9. 
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We remind the reader that the RHSs of ( l60ll(l6T] ) and thus also of (I72l)(l73l ) are invariant 
under the source scalarization transformations shown in [2] and [3l meaning that V_l x E± 
and V_L X if^and also 9 are invariant under these transformations. Therefore the 
scalar potentials \E', 9 can only depend on the unique scalarized current densities J'c and 
J'*. Whereas using the vector potential method in §[3] or the Hertz vector method in §[1] 
scalarization could only be obtained with some effort by applying appropriate source and 
gauge transformations, the scalar Hertz potentials and 9 are the natural potentials for 
obtaining scalarization since they are invariant under the required transformations. The 
final equations can then also be obtained immediately by making the RHSs of (l68l ) - ([7TI) 
zero and replacing J^, J* in (rf2l) (rf3l) by the scalarized versions J'c, J'l leading to 

Cs{e,fi){^) = -J'c (74) 
/:,(/i,6)(-9) = -J'l (75) 



However note that Ec, He and $, U are not invariant. From (17011(17111 we find 

dc 
dc J 



$ = eT W ( -u + ^ 1 rft (76) 



n = -/.I^ I [n* + ^]dt (77) 



and subsequently from (l68ll(l69ll 



9$ 

. -'H-gj-.x^ (78) 

The total fields can then be written as 

— d"^ 

E = -{v* + j^^—)c + V±e xc + V^ (80) 

— (99 

H = (v + e^ — )c + V±'^ xc + VU (81) 
at 

These expressions confirm the field transformations in [2] and [3] and they correspond term 
for term with the expressions in (I42ll(l43ll with the correspondence 

m = — = -Q = ^= 121 82 

fi± at n±_ e± at e± 

—dt n = -0* - / — dt (83) 

e± J /U± 



which was already (partially) apparent from ([TUKTfSll and (l39ll (llOll . 
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6 Gyrotropic media 

For a gyrotropic medium the transversal dielectric tensor is given by [H] 

= ex/ + jelcx/ (84) 



e± 



with a similar expression for the permittivity tensor. 

Since the scalarization of the solenoidal parts of J± and does not involve these transver- 
sal constitutive tensors no changes are needed here. However for the irrotational parts 
the reasoning leading to ( I24l)(l25l ) and ( l26l l must be extended with additional terms. In- 
troducing the gauge function Gc in the "electric" domain and as before / in the 
"magnetic" domain (fTTll we obtain the transformed stream potentials 

p = V± j udt + e-Vl + VxGc X c (85) 

m* = -^c (86) 



If e'j^ 7^ then the 2nd term on the RHS of (l85|l contains an extra term which can be 
compensated by the last term if 

a = je^/ (87) 

In this way the extra transversal term in the polarization is transformed into a longitudinal 
magnetization, which can be represented by a longitudinal magnetic current density 

'^c = -/^//-^ (88) 
Since (I2ll)(l25l) and (l26l) remain valid it suffices thus to add a magnetic current density 

/:^.^^.//| (89) 

and due to the electric/magnetic switch (m m*) in (l86l) we must also add a matching 
field transformation 



OGr 



.e 



B' - B = /^//^c = -jj--fi//uc (90) 

Again Jc + remains invariant for this additional transformation and therefore the 
formulation using the scalar Hertz potentials should still automatically lead to scalarized 
equations as has been shown in [9]. 
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Regarding the formulation using Hertz vectors extra (transversal) terms in je'j_ will also 
occur in ( ISTI ). However the functions ipe/m in (I47l)(l48ll are only part of the gauge trans- 
formation and more in general we can also add terms in G respectively in L as in (JT])® 
and (HlDdlHl) to the RHS of ((471) ((481). With in particular G = G^c (and L = L^c) we then 
obtain instead of (ISTD condition 



c X ^ ( ^ V , m 



d_ 

dc 



\l^±\ 



+je'j_c X V±^lJe + VGc x c (91) 



where = det fi± = — /i'^. This condition is fulfilled by choosing 



J 



-id_ f JJ^ 
~ dc \ 

dc e_L dc\ l/i 



/i_L 



(92) 
(93) 



The longitudinal equation (ISS]) is replaced by 



-//- 



-II 



dc 



dt 



(94) 



where —Lc is given by a similar expression as in fl93l) . We notice that except for the 
replacement of fxj^ by /i±/ the main change is the occurence of a cross-coupling 
term due to the additional gauge functions Gc, Lc which indeed mix between electric and 
magnetic stream functions. If we replace ( l4Tll by the more general expression 



r f \ V72 , ^ -1 ^ 
£.(e,;.)=V, + e//-e,--e//--- 



(95) 



then the scalarized equations for a gyrotropic medium (omitting the accents) are given 
by 



£s(e,/i)Se + je// 



d e'_L ^ /i'_L d 



dc e_L /i_L dc 
d fi'± ^ e'_L d 



dc ix± e± dc 



~dr 

dK 
dt 



Pc 



l^ll 



(96) 
(97) 



where Eg = j^IIec and = -jl-liinmc and as for the uniaxial case these Hertz vector 
components are equivalent with the scalar Hertz potentials with \1' = dJ^^/dt and —9 = 
dJ:jdt. 
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Finally we consider the formulation using ordinary but standard vector potentials. Com- 
paring (l3Ti)(l32ll with ( I47ll(l48ll the close resemblance between both formulations is apparent 
where in particular —d(j)'^*'>/dt ^ ipe/m and therefore scalarization of ( l3Tll(i32l) might also 
be possible. However whereas for a uniaxial medium the functions 0^*-' or i^e/m are suf- 
ficient for obtaining scalar equations, for a gyrotropic medium the additional freedom 
offered by the gauge function Gc in (lOTll (and L^) is needed and these functions do not 
occur in ( l3Tll (l32l) . With some hindsight we realize that transformations of the "electric" 
and "magnetic" stream potentials according to 

e"^ -p' = e'^ -p-I (98) 

e-^ -p'* = e-^-p* + S (99) 

m' = m + n (100) 

m'* = m* -n (101) 

correspond with the following transformations of the current densities 



j' = j + VxH-e-— (102) 

ot 

J =j_^.__Vx^ (103) 

These transformations are thus legitimate provided the fields are transformed according 

to ______ 

e' = E + £ H' = H + n (104) 

With the additional transformation (|102pp03p of the current densities the formulations 
based on the standard vector potentials on one hand and on the Hertz vectors on the 
other hand become fully equivalent. The former can thus also be scalarized with a proper 
choice of the derivatives of the potentials d(j)^*ydt and of the longitudinal fields 8c and 
Tic. As a result the scalar equations (I96l) (l97l) also hold for the vector potentials with the 
proper substitutions A^*^ <;=^ n(e/m)c, Pc Jc and fi//m* -v^ J* (omitting the accents). 
Also for this transformation V_l x E_i and x H± are invariant, explaining why it is 
also automatically included using the scalar Hertz potentials. 



7 Conclusions 



We have compared three "potential" methods for solving Maxwell's equations with arbi- 
trary sources in a lineair stratified gyrotropic medium where the longitudinal symmetry 
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axis (sometimes referred to as the distinguished axis) is perpendicular to the strata. In 
particular we studied the reduction of Maxwell's equations to two scalar equations (scalar- 
ization). A prerequisite for scalarization to be possible and which we have accepted 
without proof, is that the constitutive tensors should not introduce mixing between lon- 
gitudinal and transversal field components. A second condition which was often required 
in the derivations is that the transversal constitutive tensors should not depend on the 
transverse coordinates. It is perhaps interesting to note that there are no restrictions on 
the position dependence of the longitudinal properties e// and fi//. 

Introducing the conventional vector/scalar potentials or the Hertz vectors two uncoupled 
vector equations are obtained. We have shown that with the limitations already mentioned 
these equations can always be scalarized. First the sources must be scalarized: using the 
equivalence between electric and magnetic sources and gauge transformations for the 
polarizations (stream potentials) the current densities can always be replaced by unique 
current densities along the distinguished axis. The sources for the Hertz vector equations 
(the stream potentials) can be scalarized with the gauge transformations only but it is 
easier to use (l56l) and the scalarized current densities. At this stage the problem cannot 
yet be reduced to two scalar differential equations because there is still cross-coupling 
between the longitudinal and transversal components of the vector potentials or Hertz 
vectors. However using gauge transformations (see (l36ll . (l5Ti) and (l9TD ) these cross-coupling 
terms can always be eliminated and two scalar equations are obtained. When using the 
vector potentials then, for a gyrotropic medium, an additional transformation between 
electric/magnetic sources must be performed (see (I98ll - (ll04p ) and for such a gyrotropic 
medium the final scalar equations are also coupled. The two methods are found to be 
fully equivalent and the "vector potential" quantities are merely the time derivatives of 
the corresponding "Hertz vector" quantities. 

We noticed that the appropriate electric/magnetic source transformations leave the lon- 
gitudinal "total" current densities V_l x E± and V_l x H± invariant, and this obviously 
holds also for the gauge transformations. In a third method, 4 scalar Hertz potentials 
are defined based on Helmholtz decompositions (i58l)(l59l) in the transversal plane and it 
follows then immediately that 2 of these scalar Hertz potentials, \l/ and 6, are also in- 
variant under all transformations needed to obtain scalarization. These are the natural 
potentials for obtaining scalarization and this invariance explains why, when using these 
scalar Hertz potentials, the scalarized current densities emerge effortlessly. These scalar 
Hertz potentials are also equivalent with the longitudinal components of the scalarized 
traditional Hertz vectors used in the 2nd method and eventually we conclude that the 
three methods are fully equivalent. It remains to be investigated whether this conclusion 
still holds for the more complex media for which scalarization has been obtained using 
the scalar Hertz potentials [TO] . 
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